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A  bst  i  act 


Calibration  p  ocedures  for  evaluation  of  the 

acoustic  performance  of  the  Advanced  Detection  A  t  ray  (A  DA  ) 

•  v  •'  issed.  rhe  calibration  data  '  . 

1978  operation  include  measurements  of  location  and  wave¬ 
front  errors  across  the  array,  large  and  small  signal  arrav  gain, 
and  effective  spatial  signal  processing  gain  at  vt  ral  f reqt  v  , 


/> 
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i. 


IN  (  •■OIT'UC  riQN 


J  . 

Thi?  ir  orersent  of  array  piiti  is  of  coeirse  essential  to 
the  interpretation  of  the  results  of  detection  experiments 
and  background  statistics  mea  s  urements.  I  h  e  i  ain  erf'  =sis  of 
this  report  is  on  the  miesurei  ants  and  analyses  used  to 
arrive  both  at  the  array  gain,  and  at  qusnlative  estimates  of 
the  sources  cf  discrepancies  between  expected  and  measured 
per  forir.ance  Illustrative  date.  [tom  the  June.  19 ’3  seatrip 
art  presented 

The  detailed  i  nves  t  i  gat.  i  on  of  sensitivity.  frequency 
1 i ! ponse  and  directional  pattern  of  individual  hydrophones  is 
not  a  p.=  rt  of  normal  at-sea  operations  and  is  not  included 
here.  Such  hydrophone  calibrations  were  reported  in  r.ef  1. 
Lin  the  other  hand,  operations  do  include  e  cr.erk  cn  the 
status  of  each  hydrophone  and  the  selection  of  correctly 
up  c  latino  elements  for  use  in  besmf  onrinti.  II.  is  procedure  is 

d(  scribed  in  Section  II  The  me  a-  urn:  ent.  of  element  position 
e  v  i  or  5  and  wavefront  distortion,  which  affect-  the  proper 
phasing  of  the  element  outputs,  is  discussed  in  Section  111 

The  array  gain  measurement  i«  done  in  two  parts:  the 
peak  response  to  large  s  i  q  nr.  1  1  is  a  direct  measure  of  the 
ability  of  the  array  and  teamformer  to  operate  as  a  spatial 


f i Iter, 


and  is 


considered  in  S r  cion  i V.  while  the  on  —  =  x i s 
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IJ.  ELErVENT  CAL  1  BRAT  I ON  A'JD  M  l  t-CTION 

The  first  step  in  the  opn  ation  of  f DA.  once  it  is 
submersed.  is  to  reasure  the  signal  response  end  the  noise 
output  of  each  array  element.  nn<J  to  identify  ife  e  1  =nts 
which  are  not  working  properly.  The  h  earn.’ or  .  er  inputs  i  or 


the  faulty  elements  are  then  inhibited  before  any  array  gain 
measurements  c-r  processing  experiments  ai  e  started.  The 
hydrophone  survey  and  selection  process  mr.y  be  repeated 
several  times  during  the  operation,  and  a  'inal  verification 


1 un  is  made  before  the  array  is  surfaced. 


'll  e  measurements 


and  the  analysis  are  aut(  iclcd  so  that  the  selection  can  be 


completed  in  a  reasonable  time. 


Typically.  the  whole 


procedure  takes  about  three  hours. 


J  Experimental  measur ement. 

The  system  provides  fur  the  monitoring  and  digital 
recording  of  the  outputs  of  any  1  wo  selected  hydrophones.  at 
either  the  input  or  the  output  of  the  digital  filter.  Tor 
the  element  selection  process.  one  of  these  cal/mon 
(cal  i  brat  e /monitor )  channels  is  held  find  or.  a  manually 
selected  reference  hydrophone.  the  second  is  connected  to 
each  of  the  7  Z'O  hydrophones  1 1.  tuin.  Digital  filter  outputs 


ii  v  e  used 

For  the  signal  response,  the-  array  is 


oriented  t  dust d 
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U!<b  and  a  sound  source  su;  pended  from  OUH  at  approximately 
the  same  depth  as  ADA.  The  computer  generates  100  msec  tone 
bursts  at  either  1250  or  2500  Hz;  these  are  amplified  and 
drive  the  suspended  source.  '(he  cfigital  filter  is  programmed 
for  either  the  1150  to  1350  Hi  or  the  2700  to  2600  Hz  band, 
as  apprc'priate. 

The  selection  process  depends  only  on  the  response  of 
each  element  relative  to  the  selected  reference  phone*  and 
no t  on  absolute  sensitivity,  so  the  drive  level  and  signal 
conditioner  gain  are  adjusted  to  give  a  pood  s i gna 1 -t o -no i s e 
ratio,  while  avoiding  overloads  in  either  the  driver  or  the 
ax  )  ay  electronics. 

The  pulse  repetition  rate  is  varied,  depending  on  water 

✓ 

depth.  array  depth.  and  range  from  the  source  to  ADA.  to 
avoid  interference  from  reflected  arrivals.  between  pulses, 
the  computer  changes  selection  of  the  test  hydrophone, 
typically  a  1  0  to  1.  5  second  interval  is  used,  so  that  the 
scan  of  the  720  elements  may  require  twelve  to  eighteen 
minutes.  After  a  delay  to  allow  for  the  pulse  travel  time, 
cal/mon  outputs  of  both  the  test  and  reference  channels  are 
recorded.  The  record  is  made  long  enough  to  be  sure  to 
include  the  entire  pulse  arrival,  allowing  for  some  variation 
in  the  array  position  while  the  elements  art  scanned 


i 

i 


/ 


5 
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fc'l  emc-nt  noise  outputs  me  uipssured  using  the  same 
program*  but  with  the  driver  turned  off.  The  800-3200  Hi 
digital  filter  is  selected.  The  array  is  directed  away  from 
ClfVBj  to  eliminate  machinery  noise.  ’the  record  length  for 
each  hydrophone  is  0.5  seconds*  so  that  the  scan  tal.  es  six 
minutes. 

«*.  Data  reduction 

The  recorded  data  are  processed  by  programs  which  read 
the  tapes  and  produce  disk  files  containing  relative  signal 
response  or  noise  power  for  each  hydrophone.  The  signal 
response  program  uses  a  cross  correlation  filter  to  extract 
Ihe  signal  at  the  test  frequency*  in  25  msec  intervals.  The 
cm  relation  sums  are  stored  on  an  intermediate  disk  file* 
which  may  be  used  to  produce  a  calibration  report*  if 
desired  For  the  selection  procedure*  this  File  is  read  by  a 
pi ogram  which  tests  the  reference  channel  power  for 
successive  25  msec  intervals  until  the  leading  edge  of  the 


pulse  is  found  Then 

the  data  for 

a  75 

ms.ee 

i nterva 1  are 

coherently 

c  omb ined. 

and  the 

power 

1  s 

calculated  and 

converted  to 

decibels 

for  both  the 

test 

rid 

the  reference 

channels  The  difference  is  stored  in  the  relative  signal 
response  file  for  input  to  the  selection  program. 


For  noise,  the  time  samples  for  each  channel  are  squared 
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and  averaged  over 

a 

300  msec 

interval.  c :  i  d  s 

tored  in  an 

intermediate  file. 

Th  i  s 

data  is 

then  converted 

to 

noise 

power  relative  to 

th  e 

reference  phone.  The 

noise 

rr«ay 

optionally  be  normalized 

by  the  s 

ignal  response  to 

yield 

the 

equivalent  input  noise.  Again,  the  result  in  dB  is  stored  in 
a  disk  file  which  will  serve  as  input  to  the  selection 
program.  For  this  300  msec  integration  time,  the  expected 
error  in  the  noise  power  estimate  for  the  800-3200  Hz  band  is 
O.  23  dB 

8.  Selection  procedure. 

The  processed  data  files  Prom  these  runs  are  read  by  a 
program  which  prints.  or  displays  fur  the  operator.  a 
histogram  of  the  realtive  response  or  noise  measurements. 
From  these.  upper  and  lower  acceptance  limits  for  good 
hydrophones  may  be  selected.  Ihe  selection  program  generates 
a  file  which  has  three  status  indicators  For  each  hydrophone: 
one  for  signal  response  at  each  Frequency,  and  one  for  noise. 
Initially.  the  indicators  arc-  all  cleared.  As  each  of  the 
processed  data  files  is  read,  the  appropriate  indicator  is 
sc  t  for  every  element  whose  response  falls  within  the 
specified  limits.  The  final  selection  of  hydrophones  to  be 


used  by  the  beamformer  require  s  all  ttree  status  indicators 
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In  order  riot  to  reject  good  hydrophones  because  of  a 
single  bad  reading.  two  signal  response  runs  at  each 
frequency  and  two  noise  power  runs  are  made.  The  selection 
status  indicator  will  be  set-  if  either  of  the  two  readings 
meets  the  selection  criterion.  The  sssufi.ption  is  that  the 
probability  of  the  same  good  element  being  rejected  by  two 
independent  ir.eesur ements  is  small,  and  that  the  chance  that  a 
bad  hydrophone  will  be  accepted  on  at  least  one  run  on  each 
type  of  measurefi'ent  is  also  smell. 

4 .  Results  and  analysis. 

A  page  of  the  relative  response  summary  for  the  2  June 
1778  hydrophone  selection  is  reproduced  in  Rig.  II.  1.  This 
rf  port  is  not  used  in  the  selection  process,  but  is  era mined 
to  ensure  that  the  results  of  the  measurements  and  the  first 
stages  of  processing  are  reasonable.  Table  II.  I  shows  the 
response  distributions  from  which  the  operator  chooses 


ac  c  ep  tanc  e 

1  imi  ts. 

The  v  e  r  t  i 

cal  bars  indicates  the 

limits 

bit.  i  c  h 

were 

used  in  this  case. 

Two-dimensional 

h  i  s to g r 

ms  of  the  seme  runs  are  s 

r.own  in 

Rigs. 

11.2. 

These 

show  a 

r  ft  !  oriablc  consistency 

between 

T  UTl  5  . 

i .  e . 

.  most  of 

the  me 

esurement  pairs  are  clustered 

along 

the 

d  i  agonal 

.  where 

the  two  measurements  are 

nearly 

equal 

We 

can  use 

these 

histograms  to  justify. 

in  an 
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Figure  II.  1  Relative  Response  Summary 


RELATIVE  SENSITIVITY  1250  Hz 


RELATIVE  NOISE 


Figure  11.2(c)  Two-dimensional  histogram 
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informal  icay.  our  earlier  claim  that  the  relection  is  not 
likely  to  reject  a  good  hydro phone.  For  example.  at  1250  Hz. 
672  hydrophones  are  within  the  acceptance  limits  for  one  or 
both  runs,  while  ten  are  in  limits  for  Run  1  only.  and  nine 
for  Run  2  only.  This  means  that  the  probability  of  rejecting 
a  good  phone  on  a  single  run  is  roughly  10/6 72  or  0.  015.  and 
that  of  of  rejecting  a  good  phone  on  two  independent  runs  is 
only  about  0.00022. 

Figure  11.3  is  a  printout  showing  which  hydrophones  were 
inhibited  by  the  selection  procedure  for  this  data.  The 
asterisks  indicate  the  rejected  elements. 
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Figure  11.3  Killer  data  file  listing 
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111  ELEMENT  POSITION  ERRORS  AND  DOME  DISTORTION 

Errors  in  arrival  times  at  the  various  array  elements, 
relative  to  the  times  which  would  be  expected  for  a  plane 
wave,  arise  for  two  reasons:  the  true  element  positions  are 
not  exactly  known.  and  the  wavefront  is  distorted  by  the 
water  mass  enclosed  by  the  dome.  The  position  errors  are 
primarily  due  to  the  compliant  mounts,  which  do  not  stand 
exactly  perpendicular  to  the  deck;  additional  errors  may 
arise  from  varying  m.ount  lengths  and  possible  bending  or 
buckling  of  the  deck  when  the  away  is  at  depth.  The  water 
in  the  dome  is  less  saline  than  seawater,  and  does  not  come 
tc,  thermal  equilibrium  until  the  array  has  been  submierged  for 
a  long  time. 

1.  Ex  p  tr  i  menta  1  measurement. 

The  arrival  times  are  measured  using  an  e x p er i menta 1 
setup  similar  to  that  for  the  element  signal  response 
measurements.  The  primary  difference  is  in  the  element 
selection  sequence.  Again,  one  cal/mon  channel  is  held  fixed 
on  a  selected  reference  phono  while  the  other  is  connected  to 
each  element  in  turn.  The  analysis  program  will  use  the 
difference  between  the  arrival  tidies  at  the  test  and 
reference  elenients  as  a  measure  of  element  position  error 
and/or  wavefront  distortion. 
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The  ct  rot's  of  primary  interest  in  connection  with  the 
investigation  of  array  gain  arc  those  p er p end i cul ar  to  the 
plane  of  the  array.  For  these*  the  mea s ur pment  is  made  with 
the  array  oriented  toward  the  source.  However*  since  the 
array  is  not  entirely  stationary  during  the  measurements*  it 
is  necessary  to  have  a  concurrent  measurement  of  the  angle  of 
incidence  of  the  wavefront  with  respect  to  the  plane  of  the 
array.  For  this  purpose,  two  pairs  of  additional  reference 
elements  are  chosen:  one  with  elements  at  each  end  of  the 
array  near  the  horizontal  centerline*  the  other  with  elements 
at  the  top  and  bottom  near  the  vertical  centerline.  The 
selection  sequence  is  interrupted  at  regular  intervals,  and 
records  for  the  horizontal  and  vertical  pairs  are  included. 

The  signal  frequency  used  is  2700  Hz.  This  has  been 

selected  empirically  to  give  a  reproducible  maveform  with  a 
well-defined  leading  edge.  The  one  octave  digital  filter 
band  from  1600  to  3200  Hz  is  used*  since  a  narrbw  band  filter 
would  alter  the  shape  of  the  leading  edge. 
t’.  Arrival  time  determination. 

The  data  tapes  are  read  by  a  program  which  measures  the 
arrival  time  on  each  cal/mon  channel  and  stores  these  in  a 
disk  file  for  further  analysis.  The  operator  monitors  this 
}  initial  analysis  with  a  display  as  shown  in  Fig.  Ill.  1. 
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P  i  p  h  t  traces  are  shown.  four  for  each  of-  tlie  two  cal/mon 
channels.  The  upper  trace  of  each  set  of  four  is  a  51.  2  msec 
<512  sample)  "window"  of  the  data  from  the  tape.  The 
recorded  data  includes  from  one  to  four  records  of  1024 
samples,  and  the  operator  must  manually  position  the  data 
window  within  this  larger  set  of  data.  This  means  that  the 
pvopram  does  not  need  to  deal  with  multiple  arrivals.  and 
helps  the  operator  feel  needed.  Of  course.  the  pulse 
repetition  rate  must  he  set  when  the  data  is  recorded  to  keep 
surface  or  bottom  bounce  arrivals  from  interfering  with  the 
direct  arrival. 

The  second  trace  is  the  envelope  of  the  signal,  obtained 
by  e  j  p  or.ent  i  a  1  1  y  averaging  the  absolute  value  of  the  recorded 
data.  The  averaging  time  was  experimentally  adjusted  to  give 
a  reasonable  compromise  between  leading  edge  definition  and 
noise  rejection  The  envelope  data  is  scanned  for  the  peak 
level.  then  again  for  the  position  of  the  first  sample  which 
exceeds  half  the  peak.  This  gives  an  approximate  measure  of 
the  position  of  the  Jeading  edge.  The  third  trace  is  a  step 
function  whose  height  is  the  peak  envelope  amplitude  and 
whose  transition  occurs  at  the  leading  edge  of  the  pulse,  as 
determined  by  the  program. 

To  get  a  better  measure  of  the  position  of  the  leading 
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edge  of  the  pulse*  the  program  uses  an  interpolation 
algorithm  to  expand  the  time  scale  in  the  region  near  the 
edge  by  a  factor  of  eight.  Sixty -four  points  <16  before  and 
48  after  the  edge)  are  expanded  to  512  points  and  displayed 
in  the  fourth  trace.  The  effective  sample  rate  for  the 
interpolated  data  is  80  k l<2  *  which  gives  a  position 
resolution  of  1.8  cm. 

/ 

The  program  must  now  search  for  a  point  in  the  waveform 
which  may  be  reliably  identified.  for  this  particular 
waveform*  the  first  negative  peak  seems  to  be  a  good  choice. 

'I  o  make  the  measurement  reasonably  independent  of  element 
sensitivity*  the  peak  is  identified  by  comparing  the  data  to 
a  threshold  which  is  one  third  of  the  peak  envelope  level. 
This  threshold  was  empirically  chosen  to  give  as  reliable  an 
identification  as  possible. 


The  leading  edge  of  the  pulse  is  not  a  particularly  good 
point  for  determining  the  arrival  time.  More  reproducible 
results  are  obtained  by  counting  several  cycles  into  the 
pulse  to  get  beyond  the  transient  response*  which  is  not 
closely  matched  from  element  to  element.  The  program*  finds 
the  positive-going  axis  crossing  following  the  third  negative 
peak.  This  point  is  identified  on  the  display  by  a  vertical 
offset*  and  its  position  is  the  final  measure  of  arrival 
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time.  The  axis  crossing  is  used  because  the  time 
determination  is  most  precise  where  the  slope  is  a  maximum. 

3.  Time  difference  error  d  e  t  ermii  na  t  i  on. 


The 

next  step  in 

the 

processing  is 

to  compare 

the 

measured 

arrival  time 

3  "t 

e  a  c  h  element 

(relative  to 

the 

re  Terence 

hydrophone), 

to  the 

expected  time 

computed  from 

the 

nominal  position  of  the  element.  The  expected  arrival  time 
is  of  course  dependent  on  the  'direction  of  the  source 
relative  to  the  array.  so  two  passes  are  made  through  the 
travel  time  data  file.  First,  the  program  looks  only  at  the 
data  for  the  horiiontal  and  vertical  reference  pairs.  From 
these,  it  calculates  and  stores  the  components  of  a  unit 
vector  normal  to  the  wavefront.  On  the  second  pass,  the 
relative  arrival  time  for  each  element  is  compared  to  an 
expected  value  obtained  tiy  interpolating  between  the  stored 
victor  components. 

A  Results  and  analysis. 

Figure  III. 2  is  a  lifting  of  the  position  error 
estimates  for  a  typical  run.  Each  pair  of  columns  lists 
hydrophone  numbers,  from  0  to  719.  together  with  the 
estimated  position  error.  in  centimeters.  The  asterisks 
indicate  hydrophones  for  which  the  program  could  not  make  a 
position  estimate  or  which  had  been  previously  masked  out  of 
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Figure  111. 2  Position  error  estimates 
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the  array  in  the  hydrophone  select  ion  ptoccss.  Most  of  the 
errors  are  clustered  near  21  to.  Some,  however,  are  in  the 
vicinity  of  2B  cm  <a  half  wavelength  at  2700  Hz)  or  56  cm 
(one  wavelength).  The  one  wavelength  errors  result  from  a 
failure  of  the  program  to  correctly  identify  the  first  cycle 
of  the  received  pulse.  these  points  were  not  numerous  and 
were  ignored  in  the  subsequent  statistical  analysis  of  the 
data.  On  the  other  hand,  the  half  wavelength  errors  turn  out 
to  identify  hydrophones  whose  outputs  arc  phase  reversed; 
this  was  verified  by  examining  the  recorded  waveform  for 
several  such  cases.  There  were  ten  phase  reversed 
hydrophones.  and  they  were  inhibited  during  the  remainder  of 
the  operation. 

In  order  to  investigate  the  reproducibility  of  these 
measurements,  position  error  estimates  from  three  independent- 
runs  have  been  compared.  ‘the  absolute  value  of  the 
difference  between  pairs  of  estimates  for  the  same  hydrophone 
has  been  calculated,  and  a  histogram  of  these  differences  is 
shown  in  Fig.  111.3.  Ninety  percent  of  them  are  3  cm  or 
less.  The  agreement  between  the  data  for  the  three  runs  has 
been  used  as  a  selection  criterion  for  the  statistical 
analyses  which  follow.  A  position  error  is  rejected  as 


invalid  if  it  is  greater  than  20  cm  (probably  indicating  a 


Number  of  occurrences  per  1  cm  bin 


IIP!.  U  ‘<3/78 


phase  reversal  or  a  one  oevelength  measurement  error)  or  if 
it  does  not  agree  with  at  lrest  one  of  the  other  two  runs 
within  4  cm. 

Using  this  criterion)  the  mean  position  error  has  been 
found  to  be  2  12  cm.  Since  the  positions  are  relative/  this 
merely  represents  the  distance  of  the  reference  hydrophone 
from  the  mean  array  plane  and  is  not  important.  The  two 
statistics  which  are  significant  are  the  l<M3  deviation  and 
the  mean  absolute  deviation  of  the  hydrophones  from  the  mean 
aiiay  plane.  These  are  1.87  and  1.49  cm.  respectively.  The 
mean  absolute  error  affects  the  peak  beam  response,  as 
discussed  in  Section  IV.  while  the  RMS  error  is  significant 
for  the  prediction  of  loss  of  array  gain  for  small  signals, 
which  will  be  considered  in  Section  V. 

Ule  can  consider  the  apparent  position  errors  to  arise 
fi om  two  sources;  a  random  physical  position  eTror  due  to 
imperfect  mounts  and  structure.  and  a  systematic  part 
resulting  from  the  wavefront  distortion  at  the  dome.  To 
isolate  the  dome  effect,  in  table  III. I  we  have  considered 
the  array  to  be  divided  into  eleven  horizontal  slices,  each 
60  cm  high.  The  mean  position  error  for-  each  slice  is 
listed.  Also  shown  is  the  expected  error  of  the  mean,  based 
on  the  number  of  elements  in  the  slice  and  the  1.87  cm 
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Tabic  111.  1  Vertical  dependence  of  position  error 
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j  t  .Midard  deviation  for  the  whole  array.  The  standard 
deviation  of  the  set  of  11  means  is  only  O.  56  cm»  so  the 
wavefront  distortion  does  not  contribute  a  major  portion  of 
the  total  position  error. 
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IV.  PEAK  BEAM  RESPONSE  AT  HIGH  SIGNAL  TO  NOISE  RATIO 

The  ability  of  the  array  and  beamf ormer  to  correctly 
recombine  a  coherent  waveform  is  measured  by  training  the 
array  toward  the  source  suspended  from  ORB.  transmitting 
pulses  at  a  high  enough  level  so  that  the  clipped  ,b  eamf  ormer 
input  statistics  are  dominated  by  the  signal  during  the  pulse 
arrival.  and  continuously  recording  the  short  term  average 
<S1A)  beam  scan.  A  100  msec  pulse  is  used  so  that  it  always 
encompasses  the  full  50  msec  integration  time  of  at  least  one 
S'l  A  sample. 

The  data  tape  is  read  by  a  program  which  selects  and 
prints  the  largest  beam  output,  both  as  an  absolute  level  and 
os  a  fraction  of  the  maximum  possible  output.  It  also  prints 
the  beam  azimuth  and  elevation  in  both  true  and  relative 
coordinate  systems.  This  printout  is  scanned  to  find  the 
pulse  arrivals.  and  these  readings  are  recorded  for  hand 
ana  lysis. 

The  source  level  is  set  to  produce  a  s i g na 1 -t o— no i se 
ratio  at  the  hydrophone  outputs  of  approx imatel y  +10  dB. 
Measurements  have  been  carried  out  at  1  3.  1.  75  and  2.  5  kHz. 
corresponding  to  frequencies  transmitted  by  the  towed  t-rurce 
during  the  June  operation. 

Because  each  beam  response  is  a  maximum  only  in  a  single 
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direction  and  is  s  om.ewha  t  reduced  at  angles  between  adjacent 
beam  centersi  there  is  a  problem  in  trying  to  make  precise 
measurements  of  peak  beam  response  if  the  arrival  happens  to 
be  at  such  an  intermediate  angle.  It  was  expected  that  small 
random  motions  of  the  array  would  negate  this  effect  if  the 
measurement  were  repeated  several  times*  and  the  largest 
value  used.  However,  during  the  analysis  of  the  first  trial 
run  of  this  measurement,  it  was-  found  that  this  was  not  the 
case.  particularly  with  regard  to  the  beam  elevation  angle. 
The  problem  was  solved  by  creating  a  special  beam  set. 
covering  a  smaller  total  solid  angle,  with  the  beams  tightly 
packed  near  the  source  direction. 

The  peak  beam  response  was  measured  again.  and  the 
result  was  still  much  lower  than  expected,  suggesting  that 
the  phasing  was  in  error.  To  check  this  possibility.  the 
following  test  was  devised.  The  beamformer  was  reprogrammed 
by  loading  it  with  the  coordinates  of  an  hypothetical  array 
with  elements  on  a  regular  grid.  and  the  test  signal 
generator  was  enabled,  so  that  the  signal  conditioner  inputs 
were  driven  with  a  sinusoidal  signal.  This  should  result  in 
a  response  pattern  which  has  large  grating  lobes  produced  by 
constructive  interference  at  certain  angles.  In  fact,  these 
lobes  turned  out  to  be  smaller  thon  predicted. 
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By  using  different 

sets 

of 

steered  beams.  it 

was 

possible 

to  isolate  the  problem 

to 

one  of 

the 

three 

beamformer 

sections.  Fig. 

IV.  1 

shows 

the 

peak  gr 

at  i  ng 

lobe 

level  for 

beamformer  A. 

which 

wa  s 

working  correctly 

»  and 

beamformer 

Ci  which  wa  s 

not. 

The 

frequent i es 

used 

were 

chosen  to 

be  incommensurate 

with 

the 

10  kHz  s 

ample 

rate. 

since  otherwise  the  phase 

error?. 

are 

n  o  t 

random 

and 

exact 

phasing  occurs.  The  theoretical  response  for  a  uniformly 
distributed  phase  error  and  an  absolute  value  detector  is 
shown  for  comparison. 

The  error  was  found,  after  the  operation,  to  be  caused 
by  a  bad  integrated  circuit  multiplier  which  was  not  detected 
during  the  system  checkout.  The  faulty  chip  has  been 
replaced,  and  the  checkout  program  has  been  modified. 

The  peak  beam  response  was  measured  using  beamformer  A 
on) y .  The  results  are  plotted  in  Fig  IV.  2  as  a  function  of 
relative  bearing  angle.  Note  that  the  prak  beam  response 
shows  little  or  no  dependence  on  anglr.  The  theoretical 
levels  shown  are  again  calculated  for  uniformly  distributed 
phase  errors  and  an  absolute  value  detector.  The  measured 
values  are  roughly  1  dB  below  the  theoretical  levels.  The 
on-axis  measurement  at  1300  Hz  falls  very  close  to  the 


theoretical  curve. 


However.  this  sample  occured  on  the 


BEAMFORMER  C 


Figure  IV.  1  Peak  grating  lobe  levels  for  beamformers  A  and  C 
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bi  oadside  beam<  where  the  timp  delays  generated  by  the 
beamformer  are  zero  and  the  phase  error  is  not  randomly 
distributed  over  the  full  range  of  2  u  for  this  case  the 
theoretical  peak  response  could  approach  unity. 

The  measurement  of  hydrophone  position  errors  and 
wavefront  distortion  at  the  dome  was  discussed  in  Section 
111.  The  mean  absolute  error  in  the  direction  normal  to  the 
plane  of  the  array  was  found  to  be  1.49  cm.  At  2500  Hz/  this 
is  O  0248  wavelengthi  and  should  reduce  the  peak  response  by 
only  2.  5  percent.  The  reduction  is  even  smaller  at  the  lower 
Fvequenc  i  es. 
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V.  ON-AXIS  ARRAY  GAIN  FOR  SNA!  L  SIGNALS 

The  mc-a  surement  discussed  in  this  section  was  made  using 
a  sound  source  suspended  from  a  ship  at  ranges  from  5  to  12 
kiloyards.  The  transmitted  signal  was  a  sup er p o s i t i on  of 
sinusoids  at  1.0,  1.3.  1.75,  and  2.5  kHz.  The  undetected 

outputs  of  a  single  hydrophone  and  a  beam  trained  on  the 
source  were  recorded.  The  power  spectra  of  the  recorded  data 
loc-re  then  computed  and  plotted,  as  illustrated  in  Fig.  V.  1. 
This  data  can  be  interpreted  in  two  different  ways,  which  we 
shall  discuss  in  turn. 

1.  System  sensitivity. 

Here  we  want  to  consider  the  sensitivity  in  a  way 
analogous  to  the  calibration  of  a  hydrophone.  A  fundamental 
property  of  the  DIMUS  system  is  that  the  inputs  to  the 
beamformer  are  perfectly  normal i zed.  i.  e.  .  the  clipped 
signals  have  constant  power.  Clearly,  then,  the  sensitivity 
is  not  just  a  property  of  the  system,  but  depends  on  the 
nature  of  the  acoustic  environment.  For  a  signal  which  is 
small  compared  to  the  noise  <at  the  hydrophone  outputs)  the 
system  is  nearly  linear,  so  the  beam  output  can  be  used  as  a 
measure  of  absolute  signal  level. 

First,  however,  the  beam  outputs  must  be  denormalized  by 
multiplying  them  by  the  signal  level  at  the  clipper  inputs. 
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Frequency,  kHz 

Figure  V.l  Element  and  beam  spectra 
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and  corrected  by  a  frequency  dependent  gain  factor  to  correct 
For  the  beamformer  losses,  including  clipping.  In  the  ideal 
case  this  factor  would  just  be  the  n /2  clipping  loss.  The 
hydrophone  response.  in  the  direction  of  arrival  of  the 
signal,  must  then  be  applied  to  arrive  at  the  Signal  sound 
pressure  level.  These  corrections  are  not  made  by  the  system 
hardware,  but  must  be  done  after  the  fact  if.  for  example,  we 
want  to  use  the  array  to  measure  sound  pressure  levels  in 
connection  with  a  propagation  loss  experiment. 

The  plotted  spectra  have  been  used  to  calculate  the 
frequency-dependent  b eamf ormer  gain  factor.  The  calculation 
is  summarized  in  Table  V.  I.  Ihe  First  two  columns  of  the 
table  identify  the  data  record  and  the  frequency.  LB  and  NB 
are  the  line  level  and  noise  level  for  the  beam,  and  CLB  is 
the  beam  line  level  corrected  For  the  noise.  Since  the  beam 
r.  i  gna  1 -to-no  i  se  level  is  high,  the  correction  is  small.  LE. 
NE  and  CLE  are  the  analogous  levels  for  the  element;  in  this 
case  the  corrections  are  more  significant.  The  RMS  element 
output  is  available  in  the  header  of  the  spectrum  plot.  This 
is  used  to  denormalize  the  bcamFovn.er  line  level.  and  the 
result  is  listed  in  the  column  labeled  DLB.  Finally,  the 
beamformer  gain  is  obtained  by  subtracting  the  input  (CLE) 
from  the  denormalized  output  (PIB). 
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If  the  gain  factors  are  averaged  at  each  frequency/  the 
result  is  -3.4  dB  at  1.0  kHz.  -3.2  dB  at  1.3  kHz.  -5.7  dB  at 
1.  75  kHz.  and  -6.  2  dB  at  2.  5  kHz.  While  only  one  reading  was 
obtained  at  1.0  kHz>  the  readings  at  the  other  frequencies 
for  the  same  data  record  are  above  the  averages/  and  we  may 
at  least  hope  that  the  beam  we s  well  trained  on  the  target 
end  that  the  1.0  kHz  reading  is  valid. 

For  comparison/  the  theoretical  DIMU3  clipping  loss  is 
-2  dB  An  additional  correction  should  be  made  because  the 
reference  level  for  the  plotted  beam  spectra  is  full  scale 
peak  amplitude/  and  should  really  be  f<MS  amplitude.  If  we 
assume  this  correction  is  approximately  3  dB  (which  would  be 
correct  for  a  sinusoidal  wave  form)  the  expected  gain  would  be 
-5  dB/  which  lies  within  the  range  of  measured  values.  The 
lower  gain  values  obtained  at  higher  frequencies  can  probably 
be  ascribed  to  not  having  the  audio  beam  steered  exactly  on 
the  source/  which  would  affect  the  higher  frequencies  more 
strongly  because  of  the  narvower  beam.  I-or  the  04:57:20 
record/  the  frequency  dependence  is  less  pronounced;  the  beam 
may  have  been  better  trained  on  the  target. 

The  average  gain  factors  listed  above  were  used  in  the 
analysis  of  the  propagation  loss  estimates  used  in  Ref.  5. 
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£f  it  icl  processing  gain 


. Alter nativ 

e 1 y .  the 

plotted  spectra  may  be 

used 

to 

calculate  the 

process 

i  n  g  gain, 

i  .  e .  ,  the 

gain 

in 

s i g  nal -to-noi se 

ratio. 

provided  by 

the  fc csffif ormer . 

The 

s i gnal-tc  -noi se 

ratios 

at  the  o u t  p  u t  s 

of  the  h  y  d  r 

op  h  one 

and 

the  beam,  for  each  of  the 

so li t  c  e  lines. 

avf  read 

fi  cm 

the 

plots  and  compared.  Natural  1  y ,  the  spectral  lines  include 
both  signal  and  noise,  and  this  effect  must  be  corrected  to 
get  true  s i gna 1 -t o-n o i s e  ratio.  For  example,  if  the  ratio  of 
r.  i  gna  1  -p  1  u  s-no  i  s  e  to  noise  is  1  dB»  the  signal  itself  is 
about  10  dB  below  the  noise.  In  this  case,  very  small 
variations  in  s i gna 1 -p 1 u s-no i sc  mill  produce  large  changes  in 
the  result.  Since  what  me  err  trying  to  measure  is  the  gain 
For  signals  which  are  small  at  the  hydrophone  outputs,  the 
careful  selection  of  data  is  essential. 

The  spectrum  analysis  program  corrects-  the  data  for  the 
analysis  bin  width,  so  that  we  actually  have  spectrum  level 
rather  than  line  level.  However.  this  applies  to  both  the 
element  and  the  beam  data,  and  does  not  affect  their 
di  f  fere  nee. 

The  calculated  processing  gain  for  several  data  points 
at  each  of  the  four  sourer  frequencies  is  plotted  in  Fig. 
V.  2.  The  1.6  to  3.2  kHz  digital  filter  was  in  use  when  this 
data  was  taken.  and  as  a  result  the  values  obtained  at  1.0 
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«  nd  1.3  kHz  re-present  out-of-band  processing  gains  which  are 
lower  than  would  apply  if  the  system  were  operated  with  the 
low  band  or  wideband  filter.  'Ihe  reason  is  that.  although 
the  filter  reduces  the  low  frequency  components  for  both  the 
hydrophone  and  the  beam  data,  the  beamformer  clipper  noise 
will  have  a  greater  relative  contribution  than  it  would 
otherwise,  and  will  reduce  the  measured  array  gain. 

It  should  be  noted  that  this  array  gain  is  measured  with 
respect  to  the  hydrophone  outputs  and  so  dors  not  include  the 
directivity  index  of  the  hydrophones  t h t mr e 1 ve s.  The 
effective  noise  directivity  index  of  the  receiving  elements 
was  determined  by  comparing  the  outputs  from  a  hydrophone 
operating  in  the  directional  mode  and  in  the  omni  mode  with  a 
I’OO  Hz  wide  1  kHz  center  Frequency  filter  bond  These 
measurements.  made  online  with  an  RMS  voltmeter,  indicated 
that  the  effective  directivity  index  of  the  element  was 
A.  T.  dB. 
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VI.  SUMMARY  AND  CONCLUSIONS 

Element  position  measurements  confirmed  that  the 
deviation  of  the  array  elements  from  the  wavefront  of  an 
acoustic  wave  normal  to  the  deck  had  an  overall  RMS  value  of 
1  87  cm.  which  is  a  phase  error  of  .  18  radians  at  2.  5  kHi. 
The  effect  of  these  deviations  would  be  expected  to  reduce 
the  array  gain  by  0 .  1 A  d  B . 

Peak  beam  response  at  high  sigral-to  -noire  ratio  yielded 
responses  which  fell  below  the  theoretically  expected  results 
by  about  1  dB  This  difference  indicates  that  there  may  be 
some  residual  hardware  deficiencies  because  the  precision  of 
generating  a  sum  with  a  high  signcl-to- noise  ratio  should  be 
quite  high  However,  the  excellent  agreement  with  expected 
grating  lobe  response  obtained  with  the  simulator  input  would 
indicate  that  any  deficiency  was  not  in  the  beamf ormer.  but 
in  the  hydrophone  system. 

Small  signal  beam  sensitivity  mea surements  provide  the 
system  calibration  factors  required  to  infer  acoustic  line 
levels  from  beam  output  spectra.  These  factors  wer e  the  ones 
used  to  compute  propagation  lots  in  the  detection  experiment 
of  Ref  5. 


The  effective  spatial  processing  gain  of  the  array  as 
determined  by  a  comparison  of  s i mul tane ous I y  measured  signal 
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line  levels  and  broad  band  noise  spectrum  levels  on  an 
element  and  on  a  beam  output  yields  in-band  processing  gains 
of  about  24  dB.  This  is  1  dB  lower  than  the  theoretical  gain 
of  25  dB  for  a  DIMUS  processsor.  with  the  513  elements  used 
in  the  measurement.  This  1  dB  difference  is  consistent  with 
the  results  of  the  high  si  anal  -to-noi se  peak  beam  response 
nit  osur  aments 

The  overall  spatial  processing  gain  of  the  array  is  the 
sum  of  the  hydrophone  directional  gain  and  the  array  gain. 
I'4  +  4  5  =  29  5  dB 
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